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We explore the stability of domain wall and bubble solutions in theories with compact extra 
dimensions. The energy density stored inside of the wall can destabilize the volume modulus of 
a compactification, leading to solutions containing either a timelike singularity or a region where 
space decompactihes, depending on the metric ansatz. We determine the structure of such solutions 
both analytically and using numerical simulations, and analyze how they arise in compactihcations 
of Einstein-Maxwell theory and Type IIB string theory. The existence of instabilities has important 
implications for the formation of networks of topological defects and the population of vacua during 
eternal inflation. 



I. INTRODUCTION 



The introduction of extra dimensions has long been one of the most useful ideas in the theorist's tool box, providing 
insight into important theoretical puzzles ranging from the hierarchy problem to the unification of forces and a theory 
of quantum gravity. However, it has been appreciated since the early days of general relativity that there will always 
be a proliferation of low-energy effective theories resulting from compactification of extra dimensions. This idea has 
gained particular force in the context of string theory, where extra dimensions are required for theoretical consistency, 
and the program of moduli stabilization yields a vast set of low-energy effective theories, referred to as the string 
theory landscape (or simply "the landscape"). 

Classically, vacua in the landscape are stable, and drive inflation if they have positive energy density (as is the case 
in our universe if the observed dark energy is a cosmological constant). Quantum mechanically, they can decay via 
■ the nucleation of bubbles containing some new phase, a process described in many cases by the Coleman-de Luccia 
t^- (CDL) instanton [Ji]. Because the rate of bubble formation is exponentially suppressed, there is typically less than 
one bubble formed per Hubble volume per Hubble time, and therefore the inflating vacuum is never completely eaten 
up. This phenomenon is known as eternal inflation, and is a mechanism by which many different vacuum phases can 
' De populated in different spatial regions. 
t-H ' Given that the postulated existence of compactified extra dimensions is the primary motivation for considering 
landscapes of four-dimensional vacua, it is important to determine what role extra dimensions might play in the 
dynamics of eternal inflation. In the four-dimensional effective theory, the overall volume of a compactification 
^ , appears as a dynamical scalar field that inherits dilatonic couplings to the other low-energy fields. Therefore, the 
minimal model of a landscape obtained from a theory with extra dimensions contains at least two fields: one with 
properties determining the four-dimensional vacua (for example a scalar or gauge field), and the field related to the 
jj^ ■ volume, which we will refer to as the volume modulus. 

The dilatonic coupling to the volume modulus can cause a general obstruction to finding non-singular static domain 
walls and CDL bubbles connecting four dimensional vacua. Because of the coupling, domain walls can perturb the 
volume modulus out of its vacuum, in some cases making it impossible to satisfy the appropriate boundary conditions. 
This was first observed by Cvetic and Soleng [1,0] and was applied to the landscape of Type IIB string theory by two of 
the authors [j| . More recently, Yang has studied these effects in 6D Einstein-Maxwell theory [f| . A similar obstruction 
in more general multifield models can prevent the existence of other topological defects, such as monopoles Q and 
strings @4I3]. 

The static domain wall and CDL bubble solutions result from imposing a rather restrictive ansatz on the metric 
and fields, which if relaxed, may allow for additional solutions that interpolate between different four-dimensional 
vacua. To assess this possibility, we revisit and extend our previous work [5| to domain walls and bubbles assuming 
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only planar, spherical, or hyperbolic symmetry. Using a toy model, we find that whenever static solutions do not 
exist, planar walls decay into a region where the extra dimensions decompactify. Bubble walls can also seed an 
expanding region where the extra dimensions decompactify. When the phase inside of the bubble has positive or zero 
energy density (and in some cases when the interior has negative energy density), unstable walls cause the entire 
bubble interior to decompactify unless the initial radius is larger than the true vacuum Hubble constant. In this case, 
even though the initial data contains regions of two different four-dimensional vacua, at late times the configuration 
contains only the original four-dimensional vacuum and a region of the decompactified phase. We apply some simple 
bounds to assess the stability of domain walls and bubbles in model landscapes arising from Type IIB string theory 
and D-dimensional Einstein-Maxwell theory, extending the analysis in Refs. [f| Q. 

Even when static solutions do exist, they are not always stable to perturbations. Cosmologically, this is important 
for the formation of networks of domain walls or other topological defects in the early universe. Another type of 
perturbation that is particularly relevant for eternal inflation is the collision between bubbles. We find that in many 
cases, decompactification can result to the future of bubble collisions, which has important implications for the 
observability of collisions in a realistic landscape scenario (see Ref. [l3| for a review on the observability of bubble 
collisions). 

The instability of domain walls restricts the types of four-dimensional vacua that can be populated by eternal 
inflation. Although the CDL instanton cannot describe the creation of a region of some new vacuum when the bubble 
wall is unstable, we expect other types of transitions to mediate vacuum decay. If any such process has a classical 
description at late times, then even though these transitions might initially seed a region of some other four-dimensional 
vacuum, the unstable wall will quickly eradicate it, leaving behind a bubble containing only the decompactified phase. 
This might break the landscape up into sets of disconnected islands [l4|, with the members of each island populated 
only in separate realizations of an eternally inflating universe. There may even be many eternally inflating vacua that 
cannot populate any other four-dimensional vacua in the landscape. These considerations make it harder to ignore 
the question of initial conditions, a problem that eternal inflation is often purported to eliminate. The restriction on 
purely four-dimensional transitions also motivates a more general and comprehensive study of the possible dynamics 
during the early universe in theories with extra dimensions. For example, the dynamical compactification mechanism 
of Ref. |15| can populate four-dimensional vacua from a higher dimensional space, leading to a multiverse where vacua 
of different dimensionality coexist. 

The plan of the paper is as follows. In Sec. [TIJ we introduce our two-field toy model. We review the criteria for the 
existence of static domain walls in the toy model, and perform a numerical analysis of non-static solutions in Sec. Mil 
Bubble solutions are analyzed analytically and numerically in Sec. IIVI We then briefly describe the implications of 
extra dimensions for the stability of bubble collisions in Sec. [V] The stability of domain walls and bubbles in the 
Type IIB and Einstein-Maxwell landscapes are described in Sec. E2 and we conclude in Sec. IVIIl 

II. TOY MODEL 

The compactification of extra dimensions generically gives rise to many vacuum states in the four-dimensional 
effective theory. The detailed properties of the vacua and the potential landscape that connects them depend on 
the assumptions made about the fundamental theory, and we explore several representative examples in Sec. I VII 
However, without going into the specifics of these constructions, we can define a toy model that encapsulates the most 
important generic ingredients. This will allow us to distill the relevant physics governing the stability of domain wall 
and bubble solutions. In particular, our toy model consists of: 

• A field 4>, whose potential Vq(4>) has multiple minima. These vacuum states are the analog of flux vacua in the 
Type IIB or Einstein-Maxwell landscape of Sec. I VII 

• A field X) whose potential Vi(%) has a minimum at finite Xi an d approaches zero as x -> oo. This field is the 
analog of the overall volume modulus, which determines the volume of the compactified manifold. 

• A dilatonic coupling of \ to the potential Vq . Couplings of this type generically occur in the four dimensional 
Einstein frame upon compactification. 

With this set of ingredients, the potential has the structure 



V = e-"*/ M *Vb(0) +Vi(x), 



(1) 
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FIG. 1: On the left is a plot of Vo(cj)). On the right is a plot of Vi(x) for an AdS potential with a — (dashed) and an Uplifted 
potential with b = l/4a (solid). 



where we choose the specific functional forms for Vq and V\ to be 



Vi( X ) 



_ e -2 X /M x + a£ -x/M x + be -3x/M x 



(2) 
(3) 



The potential Vo ((/>), shown in the left panel of Fig. [TJ has two degenerate minima located at cj> = ±M<^. The 
potential Vi(x) is shown in the right panel of Fig. [TJ We study two limiting forms of this potential, characterized 
by a = and b = l/4a, and distinguished by their extrema and asymptotics as x ~ * 00 ■ The family of potentials 
with a = resemble models with (possibly supersymmetric) Anti de Sitter (AdS) minima and the family of potentials 
with b = l/4a resemble models where the potential is "uplifted" (and supersymmetry is broken). We refer to these as 
"AdS" and "Uplifted" potentials respectively. The region of the potential where x ~ * 00 corresponds to the volume of 
a compactification going to infinity, and we refer to this as the decompactified phase. In physical applications of the 
toy model, we are interested in situations where the local minima are situated at large x- Consequently, we expect a 
to be rather small and b to be rather large. 

For the AdS potentials, Vi(x) has a global minimum of depth V\ = —^A/27b 2 located at Xmin — M x \og(3b/2), 
and approaches zero from below as x ~ > 00 ( see Fig. [TJ left panel). For the Uplifted potentials, Vi(x) has a local 
minimum V\ = at Xmin = M x log(l/2a), a local maximum at Xmax = ^x^ S(3/2a) of height V\ — /i^.8a 2 /27, and 
approaches zero from above as x 00 • 

The full potential is shown in Fig.[2]for an Uplifted V\. There are two full minima located at {x = Xmin, 4> = ±Af^}. 
Depending on the parameters, it is possible that the minima and maxima for x can disappear at <f> ^ Ma,, opening 
up a runaway direction towards large x (as in Fig. [2]). For example, with n = 1, this occurs for 3b < fJ^/fJ^, (with 
a = 0) and a < 3fA/fi^ (with b = l/4a). The existence of this runaway direction will have dramatic implications for 
the stability of domain walls, as we now show. 



III. STABILITY ANALYSIS OF DOMAIN WALL SOLUTIONS 



A. Analytic results 



Our goal in this section will be to find static solutions where a stable planar domain wall interpolates between 
minima of the potential Eq. [TJ Neglecting gravitational effects, and assuming planar symmetry, the fields evolve 
according to 



d: 



d 2 - — 

def) 1 z ^ dx ' 



(4) 



where z is the coordinate perpendicular to the domain wall. One class of domain walls connects minima of the 
potential Vi(Xmin) to x — ► o° with (f> fixed in vacuum. Another class of domain walls interpolates between cj> = +M,p 



FIG. 2: A sketch of the full potential in Eq. [T]for an Uplifted Vi. There are two local minima, between which a runaway 
direction to the decompactified phase at x 00 opens up. 
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FIG. 3: A portion of the domain wall solution interpolating between % = Xmin and \ 00 along one of the minima for <f>. 

and 4> — —M^ with x — Xmin on either side. For this latter class of solutions, because of the coupling between <f> and 
X in the potential Eq. Q] both fields have a non-trivial profile inside of the wall. 

Domain walls of the first class, which involve only the x sector, are easy to construct. The appropriate boundary 
conditions for the equation of motion Eq. 0]are x(~°°) = Xmin and d z x{— 00) = with <j> in vacuum. Evolving in an 
AdS potential, the field reaches x(+oo) — > 00 with non-zero dzX- 1 

For an Uplifted potential, the field reaches x(+°°) - * °o with d z x(+oo) = 0. A portion of such a solution is shown 
in Fig. [31 The field interpolates between Xmin at z — > —00 and x — > 00 at z — > +00. The tension of such walls can be 



Domain wall solutions that interpolate between the vacua for <f> have a more complicated structure. We investigate 
these solutions analytically by first finding domain walls between the vacua of Vb at fixed x = X* > an d then studying 
their back-reaction on the x sector. The static domain wall solution for fixed x — X* is given by: 



calculated exactly [l|: 




(5) 



<i>Dw{ z ) = -M^tanh 



V24 




z — z*) 



(6) 



A/0 



1 



Physical examples of these types of domain walls are often supersymmetric, and thereby stable (see e.g. [T^I V 
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FIG. 4: The Euclidean potential — Vi in the Uplifted (solid) and AdS (dashed) case. Domain wall solutions must interpolate 
to and from the Euclidean maxima of the potential. 



where z* is the position of the domain wall. The thickness of the wall is approximately 



Az = -^i—e nx */ 2Mx . 



The tension is defined by [l[ 



rZ* +d 

a 4> = 2 [V^ow.xl-ViM^.Xm^dz. 

Jz*-d 



(J) 



(8) 



To capture all of the energy density associated with the wall, we must strictly take the limit where d — > oo. In the 
following, we neglect the small contribution to the tension that generically arises from the integral over V\. This is 
consistent as long as V\ (x* ) ^ Vq , which is generally a good approximation when \* ~ Xmin (stable domain walls) 
or when x* ^ Xmax (unstable domain walls). In order to keep the dependence on x* m the tension explicit, we also 



define u 



where 



-n X '/2M xz g0 that 



0> 



= 2e~ nx " /2Mx 



u* -\-d 



V du = e - nx " /2M * 



cr , 



u* —d 



00 



4^2 



(9) 



(10) 



would be the tension of the domain wall if the potential for <f> were simply V((/>) = Vo((j)). 

We now estimate how the x _ field evolves in this thin domain wall background. For static solutions (d t x = 0), the 
problem reduces to determining the motion of a particle in the inverted potential — Vi of Fig. 21 which at 'time' z = z* 
experiences an additional force from the domain wall. Substituting the domain wall solution into the equations of 
motion we obtain 



dlx 



- n *' M *Vv[4>Dw{z,X*)] + 



dV 1 
dx ' 



(11) 



Note that even when the x field starts in its vacuum (where dVi/dx = 0), the domain wall necessarily induces a 
spatial gradient. 

When the thickness of the wall is very small compared to the other scales in the equation of motion, the domain 
wall for 4> acts as a 5-function source for the x field. The influence of the dVx/dx term in Eq. [TT]on the motion of x 
is then negligible in the vicinity of the domain wall. Thus, we estimate the discontinuity in the spatial gradient of the 
X field by integrating the equation of motion Eq. [TT] across the location of the domain wall. Doing so, and using the 
tension defined in Eq. |H1 gives 



(12) 
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This relation first appeared in Ref. [3|. 

We now return to the full equation of motion for \, including the influence of V\(x)- We would like to find a 
solution where \ = Xmin on either side of the domain wall. In the mechanical analog, this means that the particle 
must start out at the maximum of the Euclidean potential — Vy at z — > — oo and end up back at the maximum as 
z — » +00. This is illustrated in Fig. [4] Because X receives a kick at the location of the domain wall in </>, staying at 
the maximum for all z is not a possible solution. Instead, \ must gain enough velocity in the positive- x direction to 
exactly balance the gain in velocity in the negative- x direction due to the kick. A solution exists if the size of the kick 
can be tuned so that it matches twice the magnitude of the rolling velocity of the particle. In this case, the trajectory 
is time-reversed, and the particle ends up at rest at the maximum, without over- or undershooting it. 

Since energy is conserved, the rolling velocity in the mechanical analog is simply given hy d z x = 
\/2\Vi(xmm) — Vi(x*)li which we must compare at different values of z (or equivalently, different values of x* ) to the 
kick Eq. [T^l There is a solution when, for some x* 7 

2V2|^(x m in) - MX*)I = ^ e "" X72Mx ' ( 13 ) 

For the AdS potentials, it is always possible to satisfy Eq.[T3] At large \i the potential goes to zero, and the rolling 
velocity is therefore approximately given by the height of the Euclidean maximum. In the limit of large b, we must 
go out to large x* to satisfy Eq. [T31 which, using the asymptotic rolling velocity, reduces to 

3V3b M x 

This is satisfied for 



2M X , 
-loe 



li\M x 



(15) 



In physical applications, it is well known that BPS domain walls exist between supersymmetric vacua (3. [l6l - l25l ] . The 
domain walls we have constructed in our toy model are analogous these types of solutions. 

The story is very different for the Uplifted potentials. In this case, the Euclidean potential — V\ has a local minimum, 
past which the rolling velocity decreases as x ~ * 00 ■ Consequently, it is possible to find situations where the kick 
always overwhelms the motion in the potential. If this occurs, the static domain wall solution does not exist. Just 
as for the AdS potentials, a static solution only exists when Eq. [T3]can be satisfied for some x* ■ Without using the 
full potential Vi(x), there are two simple tests we can apply to check this condition. These will be convenient in our 
analysis of physical models in Sec. IVI| where the potential for x takes a complicated form. 

For the first test, we can compare the rolling velocity at the minimum of the Euclidean potential (where it is 
maximized) to the kick velocity. In this case, when the inequality 



na 
2M X 



e -n Xm ^/2M x < 2 v /2|V 1 ( Xmax )|, (16) 



is satisfied, there exists a solution where x does not roll past the position of the Euclidean minimum. Substituting 
with Xmax and Vi(xmax), our toy model admits a static domain wall when 

2 (»-3)/2 3 (n-l)/2 n (tfj^j < a (2-«)/2. (1 7) 

This translates into a lower bound on a. For example, with n = 0.5, and lAM<f, = fi x M x , the critical value of a below 
which the condition for a stable domain wall solution is no longer satisfied is a cr i t — 0.18. This lower bound on a is 
an important restriction, since we should recall that a has to be small in order for the potential to have a minimum 
at large x (as is required in physical realizations of this toy model). 

Re-arranging Eq. [TSJ and noting that a x ~ \/ \V\{xmax)\M x , we see that the condition for a solution roughly 
corresponds to requiring that 

o> < °x- ( 18 ) 



Consequently, stable domain wall solutions exist only if the tension of a domain wall, separating the two vacua 
at finite x, is smaller than the tension of a x domain wall, separating the vacua at finite x from the dccompactificd 
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AdS potential (a = 0) Stable 

Uplifted potential, n > 1 Stable 

Uplifted potential, n < 1, a > a cr it Stable 

Uplifted potential, n < 1, a < a C rit Unstable 



TABLE I: For the AdS potential, planar domain walls are always stable to decompactification. For the Uplifted potential, 
stable planar domain walls are only possible for restricted values of the parameters. Recall that small a is required for the 
potential to have minima at large values of x- 



phase. If the <p domain wall has the larger tension, it is energetically favourable to form a lower-tension x domain 
wall, i.e. for the solution to decompactify. We will see in the numerical results of the next section that this is indeed 
the case. 

We now turn to the second test for domain wall stability. If the field rolls past the Euclidean minimum, it may still 
be possible to find a solution at large x* ■ I n this limit, the potential for x can be approximated by 

Vi ~ A4 ae ~ x/Mx - ( 19 ) 

If the asymptotic falloff of the kick velocity is faster than that of the rolling velocity, then it is always possible to find 
a solution at sufficiently large x* ■ Comparing yfWf to Eq. Q21 we see that this occurs for n > 1. However, when 
n < 1, then if there is no solution at small x*> there is no solution in the asymptotic region where x* — > oo. The 
parameter bounds for planar domain walls are summarized in Table HI 

When a stable solution cannot be found, the kick causes the x ne ld to overshoot its minimum and roll off to 
X — > — oo. At large negative x, the potential can be approximated as 

Vl ~ 4 e -3x/M x _ (20) 

4a 

There is an attractor solution to the equations of motion for such a potential, given by 



X{z) = log 



(21) 



From this solution, it can be seen that x ~^ ~ o° at finite z = zq, and the potential energy diverges. 

The solutions we have constructed neglect the gravitational back-reaction of the walls. Perhaps the most dramatic 
effect of including gravity is the appearance of a planar, timelike singularity in the unstable domain wall solutions, 
occurring at the value of z = zq where the potential diverges. Even for stable walls, the spacetime is altered in 
important ways. Due to the repulsive gravitational field of domain walls, the spacetime cannot be static (26l. [27|. 
Observers on both sides of the domain wall perceive themselves to be surrounded by a spherical membrane that 
accelerates towards them, reaches a turning point, and then accelerates away (see e.g. [111). However, gravitational 
effects notwithstanding, the field dynamics we have studied above should be a good approximation to the full solution 
in the limit where a <C M| and the thickness of the wall is much smaller than Mf/cr, where M4 is the four-dimensional 
Planck mass. Our criteria for the existence of domain wall solutions are therefore widely applicable. 



B. Numerical results 



In the previous section, we identified situations where it is impossible to find a static domain wall solution that 
interpolates between <j) = ±M^. The requirement that the domain wall configurations are static is fairly restrictive, 
and it is possible that there exist stable time-dependent interpolations. The relevant equations of motion in the 
absence of gravity are given by 

<fa m 

Without the static ansatz, it is in general no longer possible to find analytic solutions, and we employ a set of 
numerical simulations. These simulations should confirm the existence of static solutions over a range of parameter 
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FIG. 5: On the left is the static domain wall configuration for an AdS potential with b = 10. The fields at the beginning 
(solid black) and end (dashed red) of the simulation without artificial friction match nearly perfectly, indicating that the 
configuration is truly static. On the right, we show the numerically generated profile for x> with b = {1.25,2.5,5, 10} from 
bottom (dot dashed) to top (solid). As b is increased, the excursion of \ from its minimum increases. The prediction for the 
value of x — \* a t the location of the domain wall (Eq. ll3[) is indicated by the diamond in each case. The asymptotic expression 
for x* given by Eg. 1151 becomes accurate for b > 10. 

space, and inform us about what happens to domain wall configurations when our analytic approach predicts that static 
solutions cannot exist. We will see that in these cases, the unstable walls seed a growing region of the decompactified 
phase. 

To find static solutions, we first choose one of our toy model potentials. We start with an initial guess for the field 
configuration <f>(z,t = 0) given by Eq. [HI and x(z,t = 0) = Xmin- Adding an artificial friction term to the equations 
of motion, this initial data will evolve to the true static solution (if it exists) after a sufficient amount of simulation 
time elapses. We then check that the solution is truly static by inputting the final data from the run with friction as 
initial data in a run without friction, and seeing if there is any time-evolution. 

Applying this procedure to the AdS potentials, we are always able to find static solutions. An example for b = 10 is 
shown in the left panel of Fig. [5j where it can be seen that the initial and final configurations in the run without the 
artificial friction term are identical. The form of the solution is very similar to what is predicted by the mechanical 
analog of the previous section: the x field rolls out some distance, and then is reversed by the kick experienced at 
the location of the wall. We can perform a more detailed check of our analytic estimates by comparing the predicted 
amplitude of x inside the wall, given by solving Eq. [13] (or Eq. [15] in the limit of large b) for x* ; to the maximum 
excursion of x seen in the simulations. In the right panel of Fig. [5] it can be seen that there is excellent agreement 
over a wide range in the parameter b. 

We now turn to the Uplifted class of potentials. In the previous section, we estimated that static solutions would 
only be possible when the inequality Eq. [T7]is satisfied. For fixed n, this yields a critical value a = a cr a, below which 
a static solution should no longer be possible. Beginning at a value of a > a cr ;t, we find the improved static solution 
by the method described above. An example of such a static solution is shown in the left panel of Fig. [B] We then 
incrementally decrease a, at each step feeding in the improved static configuration from the previous step as the initial 
condition. In the right panel of Fig. 15] we plot the static profile for x{ z )- As a is decreased, the amplitude of x* inside 
of the wall increases, as predicted by the analysis of the previous section. 

Continuing to decrease a, eventually the initial conditions from the previous step cease to evolve to a static config- 
uration. This holds independently of the step size in a. Instead, as shown in Fig. field configurations interpolating 
between the vacua of <j) cause x to roll off to x — ► +co. If we set n = .5, the analytical comparison of the kick velocity 
to the barrier height (Eq. [T7| predicts that the critical value of a is a cr it — 0.18. Performing a more careful analytic 
analysis, in which the rolling velocity is compared to the kick velocity at all x* (Eq. [TB"]) . we predict a cr ;t ~ 0.168. 
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FIG. 6: On the left, we plot the static domain wall configuration for an Uplifted potential with a = 2 and n = .5. The fields at 
the beginning (solid black) and end (dashed red) of the simulation without artificial friction match nearly perfectly, indicating 
that the configuration is truly static. On the right, we plot the static profile x( z ) — Xmin, for a — {2, 1, .5, .25}. Just as predicted 
by the analytic model, the amplitude of \ inside of the cj> wall increases as a decreases. 



Numerically, we find the critical value to be between 0.171 < a cr it < 0.177. The analytic estimates are in very good 
agreement with the numerical results, and it can be seen that the simple estimate Eq.[T7]gives a good indication of the 
value of a cr it- We have repeated this analysis for various n, with similar results. For n = .25 we predict a cr it — 0.11 
(from Eq. H"3"l) . a cr j t = 0.11 (from Eq. [T7)i . and find 0.1148 < a cr j t < 0.1189 numerically; for n = .75 we predict 
a cr it — 0.18 (from Eq. H"3"|) and a cr it = 0.22 (from Eq. HTl) . compared with 0.1895 < a cr j t < 0.1962 from the numerics. 
In summary, whenever a static solution cannot be found in the Uplifted potentials, planar domain walls interpolating 
between the vacua at <f> = ±M^ inevitably seed decompactification. 

The analytic estimates in the previous section were made for an infinitely thin wall. Numerically, we can determine 
how varying the thickness of the cj> wall at fixed tension <tq affects the profile of the x-field. To do this, it is convenient 
to define a new parameter a = fi^/M^. At fixed ao, increasing a will decrease the thickness of the wall: 

/ 2 \ 1/3 

Az ex —j— . (23) 

In the limit where a — ¥ oo, the numerics should reproduce the analytic results of the previous section. 

In Fig. [51 we show the profile of the x-field for increasing values of a. As a is increased, the \ profile becomes more 
cuspy, and asymptotically approaches a fixed maximum excursion, x = X* i inside of the <f> wall. The <f> wall acts as 
a (^-function source for the x in this limit, in accord with the cuspy profile. Note also that the x ne kl becomes 
increasingly displaced from its vacuum for thinner walls. Quantitatively, the predicted value of x* f°r an infinitely 
thin wall (obtained by comparing the rolling velocity to the size of the kick in the mechanical analogue, given by 
Eq. Q2]) agrees well with the numerics: we predict (\* — Xmm)/M x — compared with (x* — Xrmri)/M x ~ 1.15 
obtained from the numerical results. 

Even in the range of a where static solutions are possible, especially near a cr it, the ability of the configuration to 
relax into a solution that separates the vacua of <f> with x remaining finite is very much dependent on the initial 
conditions. As an example, if we were to set x everywhere equal to Xmin, the value of a above which x does not 
roll off to infinity is reduced to a ~ 0.23 from the boundary for stable domain walls at a cr i t ~ 0.17. For a more 
complicated initial configuration, where the field is sprinkled randomly over different regions of the potential, the 
stability of various domains to decompactification is undoubtedly more pronounced. This "over-shoot problem" is 
certainly relevant for cosmological mechanisms of domain wall production. As an example, one can consider the 
Kibble mechanism, where different spatial regions of an initially hot universe relax into different vacua as the universe 
cools down. Without fine-tuning the initial conditions, it is difficult to imagine relaxing spatial regions onto different 
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FIG. 7: An unstable domain wall, just below the critical value of a (a — 0.125) with n = .5. On the left is a contour plot 
of 4>(z,t) and on the right we plot \{z,i). The initial configuration at t = is set by the static configuration for a — .25. 
Starting from the location of the <f) wall, a region where the \ field rolls to x 00 develops (the black region in the right 
panel), eventually spreading over the entire simulated region. The contours in the plot of x( z ^) extend from Xmin — 1.4 out 
to x — 4-8, which is well past Xmax, and in the run-away region of the potential. 

parts of the {0, \} potential landscape (at least for potentials with a near a C Hi) without causing \ to roll off to infinity 
in the regions between vacua of <fi. 



IV. STABILITY ANALYSIS OF BUBBLES 



A. Analytical results 



If we add a symmetry breaking term to the potential in Eq. [TJ 



AV 



<5/4e- 3 */^ 



M A 



A 



(24) 



then the vacua at <j> — ±Ma, have different energy, whose sign is determined by the constant A: 

A > 1, V{±M$) > 0, 
-1 < A < 1, y(+Af ) > 0, V(-AU) < 0, 
A < -1, V(±M^) < 0. 



(25) 



A cross-section of the potential at x = Xmin is shown in Fig. [9] The energy in the decompactified phase sets the zero 
of the potential. 

When the vacua on either side of a domain wall are not degenerate, finite energy solutions correspond to closed 
membranes separating the low-energy (true vacuum) phase from a background of the high-energy (false vacuum) phase. 
The lowest e nerg y configurations have spherical symmetry (in addition, deformed bubbles become more spherical as 
they expand [23|)- We therefore assume spherical symmetry (and again, neglect gravitational effects), leading to the 
equations of motion 



dt 



2d r 



dV 2 2 2d rX dV 

— , d tX -d rX = --T-, 

d(p r ax 



(26) 



where r is the radial coordinate perpendicular to the domain wall. 

Zero-energy configurations exist when the bubbles expand from rest with constant acceleration, in which case the 
spacetime possesses SO(3,l) symmetry. The orbits of the symmetry are characterized by p 2 = r 2 — t 2 . The equations 
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FIG. 8: A plot of x( z ) f° r static domain walls with {a = .18, n — .5} and varying a — {0.19, .28, 0.43, 1., 3.14, 6.73}. Increasing 
a corresponds to the c/> wall becoming thinner while retaining a fixed tension ctq. The profile x( 2 ) approaches a fixed shape 
in the limit of an infinitely thin wall, which is well reproduced by assuming a 5-function source in the equation of motion 
for \, as in the analytic treatment of Sec. Mil The predicted maximum excursion of x (f° r the above set of parameters) is 

(X* -Xmin)/M x ~ 1.11 




FIG. 9: A cross-section of the potential Vo + AV at x — Xmm for various values of A. The horizontal lines indicate the zero 
of the potential for different values of the parameter A. 



of motion in terms of this coordinate are given by 

3_M = dV 3d £ x = d V 

1 p dtp 1 p ax 

These coordinates only cover the region of spacetime in the vicinity of the bubble wall. It is necessary to analytically 
continue across the light cone r = t (where p = 0) to describe the bubble interior. In this region, the surfaces 
of constant field trace out the spacelike hyperbolas denned by p 2 = t 2 — r 2 , corresponding to the Milne slicing of 
Minkowski space. 

To obtain a bubble of true vacuum embedded in a sea of the false vacuum, boundary conditions must be chosen 
such that the fields at p — are in the basin of attraction of the true vacuum, while the fields at p — > 00 are in the 
false vacuum. To avoid a singular term in the equations of motion Eq. \27\ the derivatives of the fields at p — must 
be zero: {d p cj)(0) = 0, <9 p x(0) = 0}. Just as for the static domain walls, there are two possible classes of solutions: one 
interpolating between Xmin and the decompactified phase and the other interpolating between the vacua at <f> ~ ±M^ 
and finite \- 
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FIG. 10: Bubbles of the decompactified phase. The bubble wall profile (right panel) for a range of potentials Vi (x) + (left 
panel). As the height of the minimum increases (red dashed to dot-dashed green to solid blue), the excursion in field space for 
X becomes more and more limited. 



We begin by describing bubbles containing the decompactified phase x ~ * 00 m a background where x — Xmin and 
<\> ~ ±M0. In this case, the bubble wall interpolates between the minimum of the potential V\ (x) and some point to the 
right of the potential maximum. A number of numerically generated wall profiles and their corresponding potentials 
are shown in Fig. 1101 The range in field space over which x can travel is limited by the approximate conservation 
of energy in the mechanical analog of a particle evolving in the potential — V\. As the relative importance of the 
symmetry breaking term increases, the local minimum rises. This decreases the available energy, and thus the x field 
has an endpoint closer to the local maximum of the potential. 

With the addition of a symmetry breaking term, the vacua at <f> ~ iM^ are no longer degenerate, and we can 
search for SO(3,l) invariant bubbles containing (j> = —M^ embedded in a background where <f> = -\-M$. Just as for 
the domain walls of the previous section, both the <f) and x fields have a non-trivial profile inside the wall. In the limit 
where S <C 1, the vacua are nearly degenerate, and we can work in what is known as the thin- wall approximation. In 
this limit, the friction term in Eq. [27] causes the fields to undershoot the false vacuum unless they loiter very near to 
the true vacuum minimum for an extended range in p. The initial radius of the bubble (set by the loitering "time" ) 
is therefore very large compared to the thickness of the wall (set by the barrier-crossing "time"). Additionally, the 
equations of motion during the traversal from the true to the false vacuum become identical to the domain wall 
equations studied in the previous section (Eq. [4]) . 

Because the equations of motion reduce to those of the domain wall in the thin- wall limit, many of the results from 
Sec. IIIII carry over (with z replaced by p). Again, we can artificially freeze x at x* ■> and then study the backreaction 
of the bubble wall on the x sector. The bubble wall profile for ip is given approximately by 



4>b(z) = Af^tanh 
where the initial radius of the bubble is given by 



Ms 



(28) 



Po = |£. (29) 

The tension is identical to that of the domain wall, Eq. [9l and the difference in vacuum energy at each minimum in 
our toy model is given approximately by 

AV = 25p 4 x e- 3x — /M x . (30) 

Accounting for the dynamics of the x field, there is one important difference between the bubble and domain wall 
solutions. Because of the friction term in the equations of motion and non-zero height of the minimum, the excursion 
of x on either side of the bubble wall is limited. This was evident in our construction of the bubbles interpolating 
between Xmin and the decompactified phase; see Fig. [TU] For bubble solutions, it is therefore no longer possible to 
tune the kick size by placing the wall at arbitrarily large x* ■ This effect is most dramatic outside of the thin wall 
approximation. 

The relevant test for a stable SO(3,l) bubble solution is to compare the rolling velocity in the mechanical analog 
to the size of the kick, just as for the domain walls. Because of the limited range in x*, the parameter n does not 
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Uplifted, asymmetric potential, a > a cr it Stable 
Uplifted, asymmetric potential, a < a cr i t Unstable 



TABLE II: The stability of bubble walls is determined by the value of the parameter a. 



determine the stability of the bubble. Instead, as we summarize in Table HJ comparing the kick size to the barrier 
height (Eq. Il€>[) is a good method for determining if a solution exists. The critical value of a — a cr j t below which 
stable solutions do not exist should be nearly identical to that of the domain wall solutions. When a solution cannot 
be found, the %-field again overshoots its minimum, rolling off to \ — > — oo in finite p. 

Including gravity, the spacetime region in the vicinity of the bubble wall is described by the metric 

ds 2 = dr 2 + a(r) 2 dijf, (31) 

where dH 2 is the metric on a hyperboloid with spacelike norm. Inside the bubble, the metric is composed of hy- 
perboloids with timelike norm. The equations of motion for the fields {4>(t), x( t )} an d scale factor a(r) are given 
by 

d dV •; a ■ dV o a 2 f j) 2 x 2 \ 

where the dots indicate derivatives with respect to r. The scale factor evolves over a compact range in r between two 
zeros, each of which correspond to null surfaces. When the derivatives of both cf> and x are zer0 at both ends of the 
range in r, the zeros of a are merely coordinate singularities [l], [33]. If not, then the friction terms in the equations 
of motion cause a curvature singularity to develop 0, ■ This will be the case when the size of the kick overwhelms 
the rolling velocity, and the x held is forced to overshoot its minimum. For such solutions, there is no region of true 
vacuum inside of the bubble. Instead, the expanding bubble wall contains an expanding spherical timelike singularity. 

The field dynamics we studied above are valid when the radius of the bubble (defined in the thin- wall approximation 
by Eq. |2T)1) is much smaller than the Hubble radius defined by the energy of the false vacuum. When gravitational 
effects arc important, there is good reason to expect that stable non-singular solutions are even harder to find. In this 
limit, the Hubble friction terms in Eq. 1321 become important, causing the field to roll in the steepest direction (this is 
similar to what happens in slow- roll inflation), which inside of the <j) wall is towards x ~^ —oo. The same is true for 
thick walls in the absence of gravity. The result is again a timelike singularity. 

B. Vacuum transitions 

Analytically continuing Eq. [3TJ the Coleman-de Luccia instanton is obtained. This solution to the Euclidean 
equations of motion (Eq. I32[) mediates the nucleation of true vacuum bubbles from a false vacuum phase [l] . The 
singular bubble solutions studied above do not analytically continue into finite action Euclidean instantons. Therefore, 
the presence of the x -ne ld can cause the rate for certain transitions out of the false vacuum to go to zero, as first 
noted in Q. The nucleation of bubbles containing the decompactified phase, which yield finite action instantons, are 
always allowed. 

The lack of certain 0(4)-invariant instantons does not preclude the existence of an instanton with less symmetry. For 
example, 0(3)-invariant instantons [3ll. |32| can create a localized region of true vacuum. However, if the underlying 
potential does not allow for a stable bubble wall, then as we show in the next section, our numerical results indicate that 
when the vacuum energy of the true vacuum is non- negative, the entire bubble interior goes over to the decompactified 
phase. The only processes that can evolve to a configuration which at late times contains multiple vacua would seem 
to be a horizon-sized fluctuation out of the false vacuum (larger than both the true and false vacuum horizon). This 
is similar to the Hawking- Moss transition [33[ or the nucleation of a false vacuum bubble [34[, which must be larger 
than the horizon size in order to overcome the encroaching true vacuum exterior. It is unclear how to interpret such 
events, or even whether they occur at all. In the absence of such fluctuations, this implies that lasting regions of some 
vacua will not be created by vacuum transitions. 

In the physical examples of Sec. lVIl we will consider analogs of the Brown-Teitelboim [35j j mechanism for membrane 
nucleation. In these models, there is a charged fundamental brane (the thin-wall limit of the solitonic <f> walls in our 
toy model) that inherits a coupling to the x-field. The stability analysis is identical to that presented above, since the 
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FIG. 11: The evolution of (j> (left panel) and \ (right panel) for a stable bubble configuration in the {r, t} plane with a potential 
specified by the parameters {A — 1.0, a = 0.19, 8 = 0.05, n = 0.5}. The panels show a cross section of the true vacuum bubble, 
which is separated from the surrounding false vacuum by an expanding domain wall. 

branes are characterized by their tension, which in turn determines the effect on the x field. Further, in a background 
field, we expect there to be some process that creates these charged objects, described by a thin-wall CDL or 0(3)- 
invariant instanton. Such a process would be analogous to the Schwinger pair production of electromagnetic charges 
in an electric field. Again, if the wall is unstable, the resulting configuration will not contain a region of true vacuum 
at late times, just as in our toy model. 

C. Numerical results 

In order to investigate spherically symmetric solutions that do not have a boost symmetry, we again employ a 
set of numerical simulations. The relevant equations of motion are given in Eq. [2H We focus on potentials with 
6 -C 1, where the thin-wall approximation holds. In this case, as discussed above, the shape of the field profile can 
be approximated by the domain wall configuration found from the potential with 5 = 0. We therefore expect to find 
similar run-away behavior for a < a cr i t , although the development of the instability will be more involved because of 
the potentially non-zero energy in the true vacuum. 

Choosing a particular potential, we first set 6 = and find the improved domain wall configuration using the 
methods described in Sec. IIIIBI We then translate this solution, such that the center of the wall is located at a 
radius given by Eq. 1291 This is used as the initial condition in a simulation with 5 -C 1, and evolves smoothly into an 
expanding bubble solution. One such stable bubble is shown in Fig. [TT] 

Because the symmetry-breaking term changes the vacuum structure of the potential slightly, the improved domain 
wall configuration we input for x does not initially interpolate exactly between the potential minima. This intro- 
duces oscillations about the false vacuum, which are apparent in the numerical plots shown in this section. These 
disturbances are a small perturbation, and will not affect the conclusions of our analysis. 

Decreasing a towards the critical value and repeating this procedure, the bubble configuration becomes unstable, 
seeding a region where the x field rolls off to the decompactificd phase. It is difficult to find the exact SO(3,l)-invariant 
profile for stable bubbles, and thus difficult to determine the sharp value of a cr ; t , although it is close to that of the 
corresponding planar domain walls. The sensitivity to small departures from the exact profile indicates that walls are 
not stable to small perturbations for potentials with a near a cr it- 

When it exists, the way that the instability develops depends on the vacuum structure of the potential. If the 
true vacuum is negative (A < 1 in our toy model), it is of lower energy than the decompactified phase at x ~~ * 00 ■ 
Therefore, if a region of the decompactified phase forms, it will naturally be repelled from the true vacuum interior. 
An example of this behavior is shown in Fig. 1121 where it can be seen that the interface between the true vacuum 
and the decompactificd phase is accelerating out of the bubble interior. Taking a spatial cross section through this 
configuration, the bubble wall contains an expanding shell of the decompactified phase surrounding a growing true 
vacuum center. 
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FIG. 12: The evolution of (j> (left panel) and \ (right panel) for a bubble configuration in the {r, t} plane with the parameters 
{A — 0.1, a — 0.18, S — 0.1, n = 0.5}. The unstable wall creates a region in which x rolls off to +oo (black), as shown in the 
right panel. Since this is a higher energy phase than the negative energy bubble interior, the true vacuum region grows, leading 
to a nested configuration with a true vacuum core (light blue), surrounded by a region where \ — >• oo (black), surrounded by 
the false vacuum (light blue). 

The dynamics of the instability can also be affected by changes in the behavior of the (j> wall, which lies in the region 
undergoing decompactification. In the vicinity of the wall, the symmetry-breaking term in the potential is becoming 
increasingly irrelevant when compared with Vo- A decrease in the vacuum energy difference across the <f> wall (as 
compared to the tension) can change the critical radius for an expanding bubble, and actually reverse the expansion 
of the wall. An example of this behavior is shown in Fig. 1131 which occurs in our toy model for small values of 5. 
Taking a spatial cross section through the simulation, the bubble wall spawns a shell of the decompactified phase, 
which then implodes, eventually eating up the entire core of true vacuum. Because the critical radius for expansion 
diverges as x -> oo, the turn-around of the 4> wau might always occur, although we have not been able to confirm 
this due to the limitations on the size of our simulation. 

Potentials with a negative energy true vacuum give rise to a big-crunch inside of the bubble [3fl l36j . When the 
wall is unstable, we saw above that the decompactified phase exists in a region separating the bubble interior from 
the wall. It would be interesting to determine how the singularity develops in the presence of an such an unstable 
wall. 

When the true vacuum has zero or positive energy, it is greater than or equal to the energy of the decompactified 
phase, and the decompactifying region accelerates into the true vacuum. The dynamics of the <p wau can a l so drive 
the growth of a region containing the decompactified phase. For a bubble near the critical radius, as in Fig. ITU 
the growing region of the decompactified phase and associated inward moving <f> wall eventually remove the entire 
region containing the true vacuum. Making the initial radius of the bubble larger, as in Fig. 1151 the <f> wall can move 
outwards (at least for some period of time), but the region containing the decompactified phase again accelerates 
inwards, removing the region of true vacuum from the spacetime. The results from a number of simulations imply 
that unstable bubbles separating the vacua at (f> ~ ±M^ always evolve to bubbles containing the decompactified 
phase. 

If we were to include gravity, the Hubble expansion of a positive energy true vacuum interior provides an extra 
restoring force against the encroaching decompactified phase. If the initial radius of the bubble were comparable to 
the true vacuum horizon size, and the instability developed on a time-scale comparable to the true vacuum Hubble 
time, then it is possible (by causality) that the true vacuum phase inside of the bubble is not completely eaten up. If 
the bubble radius were larger than the true vacuum horizon size, then the bubble interior is guaranteed to survive. 
However, since the construction of such configurations would require acausal correlations, it is not clear that they are 
physically well-defined. 
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FIG. 13: The evolution of cf> (left panel) and \ (right panel) for a bubble configuration in the {r, t} plane with the parameters 
{A — 0.1, a — 0.18, S — 0.05, n = 0.5}. In this case, the reduction in S changes the dynamics of the <^>-wall in such a way that it 
implodes into the true vacuum. The region of true vacuum is eventually removed from the spacetime. 




FIG. 14: Here, we show a contour plot of a configuration for cj> (left) and x (right) that initially contains a true vacuum bubble 
(where (f> = +M^) embedded in a surrounding false vacuum (where <j) — —M^). The parameters of the potential are given by 
{A = 1, a = 0.19, 8 = 0.05, n = 0.5}. When the phase inside the bubble has zero or positive energy, situations with an unstable 
wall lead to the entire bubble interior running off to \ — > oo (black), as shown in the right panel. In addition, for this set of 
parameters, the <fi wall implodes into the bubble. 



V. COLLIDING BUBBLES WITH DILATONIC COUPLINGS 

In the numerical examples of the previous sections, we have seen that domain walls in compactified theories can be 
relatively fragile. An additional perturbation to consider, which is very relevant in the context of eternal inflation, is 
the collision between bubbles. Although the phase transition from the false to the true vacuum does not complete if 
inflation is eternal, there are an infinite number of collisions with any given bubble. This takes on direct observational 
relevance because within the picture of eternal inflation, our observable universe could lie within one such bubble, and 
collisions might leave signatures in cosmological probes such as the CMB (for a review on the observability of bubble 
collisions, see Ref. [Hj]). 

If eternal inflation is driven by a potential landscape obtained from a theory with extra dimensions, then the volume 
modulus can play a dynamical role in the collision between two bubbles. Each bubble is "dressed" with a non- vacuum 
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FIG. 15: A contour plot of a configuration for 4> (left) and \ (right) that initially contains a true vacuum bubble (where 
4> = +M^>) embedded in a surrounding false vacuum (where (f> = —M$). The parameters of the potential are given by 
{A = 1, a = 0.22,5 = 0.05, n = 0.5}. Using an initial bubble radius much larger than the critical radius, the (f> wall can be 
made to expand. However, the region where \ ~ ► 00 (black region in right panel) still removes the entire true vacuum core 
from the bubble interior. 



configuration of the volume modulus, as we found in the toy model of Sec. IIVI These lumps of the volume modulus 
are accelerated with the bubble wall, and since the potential is non-linear, they interact upon collision (this was also 
suggested in Ref. @). Bubble collisions can cause a very large perturbation to the configuration of the volume modulus 
in the collision region, possibly leading to decompactification. This is a generalization of the "classical transition" 
mechanism of Ref. [37] to multi-field potentials. 

We assess the ability of collisions to seed decompactification by estimating the energy released when two wave- 
packets of the volume modulus collide. The energy density in each of the wave-packets will be approximately given 
by the gradient energy estimated in Sec. IIII Al When the bubbles are large, they are approximately planar at the 
location of the collision, and transforming to the frame where the wall is expanding, we can replace d z \ — > 7(i)9 z x 
where ^(t) is the Lorentz factor of the expanding wall. Using energy conservation, as in Eq. 1131 the gradient energy 
associated with each of the wave-packets of \ in the center of mass frame is of order 

7 2 (d,x) 2 ~7 2 Vi(x*). (33) 

When the wave packets attached to each wall collide, some fraction / of the energy density thermalizes due to 
the anharmonic nature of the potential. Both the importance of the anharmonic terms and the gradient energy 
increase with x* ~ Xmm- It is also possible that there is energy transfer between the </> and x fields, since they are 
coupled. If the energy density thermalized is greater than the barrier height, then the collision has the ability to cause 
decompactification. Neglecting gravity, the condition for no decompactification is 



(34> 

where / is the fraction of collision energy that thermalizes. The center of mass energy of the collision increases with 
initial separation, and in the limit where the initial size of the bubble po is much smaller than the separation s, the 
lorentz factor can be written as 7 — s/po. The condition Eq. 1341 therefore amounts to a restriction on the kinematics 

of the collision. 

Alternatively, using Eq. Eq. O noting that a x ~ \/Vi(xmax)M x , and neglecting factors of order one, we can 
re-write the condition for no decompactification in terms of the tensions of the <fi and X walls: 

7<7^- (35) 
/ 1 ^ 



In this picture, the collision can supply the energy necessary to transition from a stable <fi wall to x walls. This 
suggests, from a slightly different perspective, that decompactification will result from sufficiently energetic collisions. 
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FIG. 16: The collision between two stable bubbles with {a — 0.5, n = 0.5, A = 2, 5 = 0.02}, which does not lead to decompact- 
ification in the future of the collision. The panels show the <f> (left) and x (right) evolution of one of the bubbles, which collides 
with the periodically identified boundary of the simulation. After the collsion at t ~ 60, \ an d 4> relax to the true vacuum. 



In order to verify that collisions can indeed cause decompactification, we have run a number of numerical simulations. 
Following e.g. Ref. [l3|, we transform to a set of hyperbolic coordinates appropriate to the SO(2,l) symmetry of the 
collision spacetime. By imposing reflective boundary conditions (i.e. imposing symmetry in the fields about each 
boundary) , we model the collision between two identical bubbles by colliding an expanding bubble with the boundary 
of the simulation. The bubble starts at rest and grows from some fixed initial radius. We vary the kinematics of the 
collision by changing the spatial size of the simulation. 

Generally, for a wide range in the parameters of the potential, {a, n, A, <5}, we find that increasing the separation s 
between the two colliding bubbles causes the \ field to get a larger kick to the future of the collision. For some range 
of kinematics, the collision relaxes back to the true vacuum, and decompactification does not occur. An example is 
shown in Fig. 1161 

Increasing the separation between the colliding bubbles, we observe two different behaviors depending on the value 
of A. If A is close to one, so that the true vacuum is nearly zero energy, then as shown in Fig. [T71 a large pocket of 
the decompactified phase forms and grows, but eventually re-collapses. Increasing A, there is a separation past which 
a growing pocket of the decompactified phase is formed, which then accelerates into the bubble. An example is shown 
in Fig. [TH1 The energy of the true vacuum must be sufficiently large to guarantee that the pocket of decompactified 
phase accelerates into the bubble. This may render vacua with a small cosmological constant more stable towards 
decompactification from collisions (note that in our universe, the relevant energy scale is that of inflation, and not 
the present day cosmological constant). 

Quantitatively, we can compare our numerical results to the prediction Eq. 1341 Running a number of simulations, 
we find that Eq. [MJ accurately predicts the lorentz factor necessary to produce decompactification for an efficiency 
somewhere between 0.1 < / < 0.5. 

In flat space, the bubble separation can increase without bound, and there is always a set of kinematics that leads 
to decompactification, no matter how inefficient the thermalization process is. In a de Sitter (dS) background, the 
typical separation for collisions is of order the maximum allowed by causality [38j, s ~ Hp 1 . This sets an upper 
bound on the achievable 7 factor, of order 7 ~ (Hppo) . Thus, the ability of a collision to decompactify is limited 
by the size of Hppo, which also parametrizes the importance of gravitational effects in the formation and evolution 
of bubbles. When gravitational effects can be neglected, the collisions will in many cases be sufficiently energetic to 
cause decompactification. 

Our cosmology cannot exist to the future of collisions that cause decompactification in a region that accelerates into 
its surroundings. This simple fact has very important implications for assessing the feasibility of observing cosmic 
bubble collisions during eternal inflation. If decompactification is a fairly generic outcome of collisions, then there will 
be restrictions on the strength of possible signatures and the number of collisions we might expect to be in our past. 
Given that the primary motivation for considering landscapes of vacua that drive eternal inflation is the possible 
existence of extra dimensions, this effect is quite important, and deserves further study in the context of specific 
models. 
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FIG. 17: The collision between two stable bubbles with {a — 0.5, n = 0.5, A = 1, 8 = 0.05}, which results in decompactification 
in some finite region to the future of the collision. The panels show the <f> (left) and x (right) evolution of one of the bubbles, 
which collides with the periodically identified boundary of the simulation. After the collsion at t ~ 100, x gets kicked to the 
runaway region of Vi in a localized region. This pocket of the decompactified phase (black) expands, but ultimately re-contracts. 




FIG. 18: The collision between two stable bubbles with {a — 0.5, n = 0.5, A — 2, S = 0.02}, which results in decompactification 
in the entire future of the collision. The panels show the <j> (left) and x (right) evolution of one of the bubbles, which collides 
with the periodically identified boundary of the simulation. After the collsion at t ~ 75, x rolls off to infinity (black). 



As an interesting aside, there is one sense in which we could exist to the future of a decompactification-inducing 
bubble collision. If eternal inflation were to occur in three dimensions, then bubble collisions might seed decompacti- 
fication into a four-dimensional universe. There may be interesting observational consequences of such a scenario. 



VI. DILATONIC DOMAIN WALLS IN COMPACTIFIED THEORIES 



In this section, we describe the link between the toy models and more realistic physical models with compatificd 
dimensions. In particular, we describe how dilatonic couplings between fields arise in theories with extra dimensions. 

After deriving some general properties of these theories, we focus on compactifications of Einstein-Maxwell theory 
and type IIB string theory. The former compactifications are relatively simple and under good technical control. 
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Compactifications of string theory are more involved, and viable models with four-dimensional dS vacua have only 
been constructed recently. On the other hand, since string theory is at the moment the best candidate for a theory 
of quantum gravity, these models are of high interest. 

In both the Type IIB and Einstein-Maxwell examples, we begin with a D-dimensional action of the form 



M D ~ 2 



d D xy/-g {D *>n (D) + s fie i ds + s s 



(36) 



Here Md is the Z?-dimensional Planck mass (with 8itGd = Mp~ 2 ), the tilde denotes that each quantity is evaluated 
in the _D-dimensional Einstein frame, Sfi e id s includes the matter fields present in the theory (e.g. gauge fields and 
bulk scalars), and S sources includes any branes or other sources that may be present. The fields and sources depend 
on the theory under consideration, and we defer discussion of these pieces of the action to the particular examples in 
following sections. 

We generally assume that the Z?-dimensional space-time factorizes as .M^ = M.4 x w i} q , where Q q is a compact 
q-dimensional manifold and w indicates that warp factors might be present. The metric is block-diagonal, and in the 
absence of warping, we will begin with the ansatz: 



ds 2 = exp 



2q X(x) 



q + 2 M 4 



g$(x)dx' l dx v +M-2 i exp 



x( x ) 



q(q + 2) M 4 



g$ n dy m dy n . 



(37) 



The field x( x ) is the volume modulus of fl q . There may also be non-trivial warping in different regions of the compact 
manifold, as will arise in the Type IIB examples below. In this case there are a number of subtleties which we comment 
on in Sec. EH1 

Inserting this metric ansatz into Eq. 1361 and integrating over the compact space, we obtain 



S = 



M 2 1 



+ 



(38) 



where 



M A = M 



D 



d q yVg^> 



1/2 



M D Vol(Q< 



a/2 



(39) 



defines the four-dimensional Planck mass. The potential for x depends both on the field content of the higher dimen- 
sional theory and the assumptions made about the form of the compact manifold. We assume that the dimensional 
reduction is self-consistent - namely, that the Kaluza-Klein modes can be integrated out and the volume of the 
compact manifold is somewhat larger than M D q . 

We now consider a number of examples where the dimensional reduction can be carried out explicitly. 

A. Einstein— Maxwell theory 

Perhaps the simplest example of a theory possessing dS vacua is the Einstein-Maxwell system with a positive 
cosmological constant compactified on a g-dimensional sphere (recent papers describing this model include |6lll5l.l39j). 
The D dimensional action for the fields in this example is given by 



M 



D-2 
D 



d D xV Z ¥^ 



(40) 



where F q is a q-form magnetic flux wrapping the compact manifold (equivalently, we could add the dual 4-form electric 
field strength), and A > 0. The field equations for the g-form are satisfied for 



F q = Q sin 9 " 1 6>i . . . sin0 g _id0i . . . A d0 q . 
In the dimensionally reduced theory, these sources give rise to a potential for x given by 



(41) 



V(x) 



Mi 



2Voi(sq) 



q(q — l)e ' v ' M * + 



2A r -A(g)^ , Q_ -3A(g)^ 



(42) 
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where the first term is due to the positive curvature of the compactification manifold, we have introduced X(q) = 
( g ^?2) ^° snT1 plify the expressions, and Vol(S f<? ) is the volume of a unit g-sphere 



Vol(5 9 ) 



2 >+i)/2 



(43) 



Depending on the charge, there can exist positive, negative, or zero energy minima of this potential, giving rise to 
four-dimensional dS, AdS, or Minkowski vacua. A sketch of the potential for various Q is shown in Fig. [T9l 

We now turn to S sources . A fundamental m + 2-brane in the D-dimensional theory can give rise to a 2-brane in four 
dimensions, with to legs of the brane wrapped on a cycle of the compact manifold. For the spherical compactifications 
under consideration, stable configurations only exist for an unwrapped (m = 0) or fully wrapped [m = q) brane. The 
to + 2-brane naturally sources electric or magnetic flux, which can decompose upon dimensional reduction into legs 
living on the compact portion of the manifold, the non-compact portion of the manifold, or both. For example, a 
2-brane can source a q-form magnetic flux living on the sphere or a 4-form electric flux living in the four non-compact 
dimensions. A wrapped q + 2-brane sources both a g-form flux on the sphere and a 4-form in the non-compact 
dimensions. 

Moving across a 2-brane, the effect of the charge is to shift Q — » Q — q. Moving across a wrapped q + 2-brane, the 
effect is to shift A — > A — q. 2 This changes the vacuum structure of the potential for \, and therefore the membranes 
separate four-dimensional vacua with different vacuum energy (hence, the membranes must be spherical to have finite 
energy). In analogy with Schwinger production of electron-positron pairs in a constant electric field, spherical charged 
membranes are nucleated from the background fields [35| used to stabilize the volume modulus (see Ref. [39] for a 
recent discussion of membrane production in the Einstein-Maxwell model). 

The source action in the D-dimensional Einstein frame, including a coupling to the 3 + TO-form gauge field A, is 
given by 



~ j3+m. 



q I A, 

Is 



(44) 



where a is the tension of the brane, q is the charge, and the integrals are over the world-volume S of the brane; 7 is 
the worldsheet metric. 

The dimensional reduction to the four-dimensional Einstein frame is performed using the metric Eq. 1371 In general, 
the addition of sources deforms the spherical compactification manifold, introducing warping into the metric (see 
e.g. [40| for a discussion of warped compactifications of this type). We neglect this effect, which should be a good 
approximation in the limit where Q 3> q. This also implies that the vacua on either side of the brane are nearly 
degenerate, and so we focus on the role of the tension, which in the four-dimensional Einstein frame contributes an 
action 



a f ^d?£,d m t = Ml 

JT, 




3q - 2m x(0 



(45) 



2 A D = q + 4-flux, as the one sourced by the q + 2-brane, acts effectively as a cosmological constant in a D-dimensional theory. 
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Here 7 is the worldsheet metric defined with respect to the four-dimensional Einstein frame metric g$, Tl m is the 
period of the m-cycle that the brane wraps (for the purpose of illustration, we have allowed for any number of legs 
of the brane to be wrapped; for m — q, IT m = Vol(fi q )). The exponent in the coupling to x is negative definite since 
m < q. 

We can now make connection with the toy model described in Section |TT] (which the reader should now recognize as 
arising from six-dimensional Einstein-Maxwell theory). Due to the tension term in the action, the potential for the 
volume modulus x gets an extra contribution at the location of the brane, inducing a spatial derivative approximated 
by 



-Mi 



M 3+m 



n, 



\ 3q — 2m 



exp 



3q — 2m x* 
^/2q(q + 2) M4 



(46) 



where the prime denotes a derivative with respect to z or r, depending on the assumed symmetry of the wall. As 
described in Sections IIIII and HVl the existence of static planar or SO(3,l)-invariant solutions including the brane and 
X field depends both on the asymptotics and parameters of the potential for x- When A = 0, the local minimum is 
AdS and the potential for x approaches zero from below, resembling the AdS type potentials in our toy model. In 
this case, as we found in Sec. IIIII there is always a stable domain wall solution. For A > 0, the local minimum can be 
dS and the potential approaches zero from above, resembling the Uplifted class of models. For such potentials, there 
is only a non-singular static or SO(3,l) invariant solution when the spatial derivative induced by the brane is smaller 
than the possible gradients induced by the potential. 

Applying the first analytic estimate discussed in Sec. IIII) we compare the change in the spatial gradient of the field 
X (Eq. I46p to the height of the potential maximum for \, The location and height of the maximum are given by [l5j 



y/q(q + 2) 
2V2 



log 



M 2 



V( X n 



Mi 



2A 

Vol(S ( l)(2 + g ) 1 +?/2( g _ 1)9/2 \J^T 



1+9/2 



(47) 
(48) 



Substituting with x* = Xmax, the ratio of the kick to the rolling velocity in the potential V\ scales like 



(Ax') 



n„ 



2V2K 



Vol(S*9) Mf+ m 



2A 

Ml 



(9-™-l)/2 



(49) 



The first factors in this expression are typically order one numbers, since ITo^ = (1, Vo^S* 9 )) (unwrapped and wrapped 
cases respectively) and the tension of the branes should scale with Mp. 3 In the last factor, we expect A/Mf, to be 
small for the semi-classical approximation to be valid. Thus, the magnitude of the ratio depends on the values of q 
and m, i.e. on how the branes wrap the extra dimensions. 

When q — m — 1 > 0, the kinetic energy is always small compared to the potential energy, and domain walls are 
expected to be stable. For q — m — 1 < 0, the ratio is large and it seems difficult to obtain stable domain walls. 
Maximally wrapped branes are therefore always unstable. Indeed, the asymptotic falloff of the gradient due to the 
brane is always larger than the asymptotic gradient due to the potential for x, ensuring that stable solutions are 
truly impossible to construct. Exceptions to this rule are only possible if the tension is tuned to be sufficiently small 
compared with which seems difficult to obtain for gravitating branes that can be described by the Einstein- 
Maxwell model. 

From this example, we infer that when a brane wraps some internal cycle, it is more difficult to find stable domain 
walls. This is due to the fact that the brane tension in the four-dimensional Einstein frame falls off more slowly 
with increasing x if the brane has one or more legs in the compactification manifold. In a model where the domain 
walls are unstable, their nucleation cannot proceed via the standard 0(4)-invariant instanton. However, one might 
expect there to be some process by which membrane production can still occur via some 0(3)-invariant [3l|, HI] or 
constrained instanton [4lj |. In this case, once formed, if the interior of the membrane has zero or positive energy density 
it decompactifies as described in Sec. IIV CI Therefore, some vacua cannot be populated by membrane nucleation in 
the Einstein-Maxwell model. 



3 If these are gravity solutions, then the horizon associated with the branes scales with the tension 
be under good theoretical control, and we therefore expect that <r > M^ +m . 



sub-planckian horizons would not 
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B. Compactifications of type II string theory 



There exist a multitude of compactifications of string theory to four dimensions. M theory, heterotic and type II 
string theories can all result in effectively four-dimensional models. 4 However, matching these models with phenomeno- 
logical requirements is non-trivial. In particular, obtaining four-dimensional dS vacua has proven to be difficult. In 
this section, we focus on compactifications of type IIB string theory, where currently the most reliable constructions 
of dS vacua have been performed. 5 We only recall the main features of these models (readers interested in the details 
and a more extensive list of references are referred to one of the many reviews on the subject, e.g. (42I . IHol . [5ll|). 

Type II string theories contain a wide variety of stable branes, sourcing different higher-dimensional fluxes. Allowing 
these fluxes and branes to pierce and wrap cycles in a compact manifold can stabilize the shape and size of this 
manifold. For type IIB compactifications we can choose the internal manifold to be conformally Calabi-Yau, in which 
case it is described by complex structure and Kahler moduli. The volume modulus x defined in Eq. [37] is a Kahler 
modulus. The complex structure moduli are analogous to the (j> field in the toy model of Sec. HU 

As was mentioned in the previous section, the addition of branes and fluxes generally leads to non-trivial warping 
of the com pac tification manifold. Warping is of phenomenological interest for generating hierarchies in particle 
physics [52 55], and is a key ingredient required to stabilize moduli in some scenarios (5o| . There are a number of 
subtleties that come with identifying moduli and determining their dynamics in a theory with warping [57l - l60j . We 
are mainly interested in the properties of complex structure moduli and the overall volume modulus, and so it is 
sufficient to use the metric ansatz of Rcf. 1601: 



ds 2 = 



,2Q(c) 



[e-*Mv) + c(x)]' 
where c is the overall volume, related to \ by 



w (g${x)dx»dx» + 2d ] Bd^df) + [e- iA °^ + c(x)] 1/2 g^(y)dy i dy^ (50) 



c(x) = exp 



2 x(s) 
3 M 4 



(51) 



and Aq (y) is a warp factor associated with the metric g^ . B is a compensator field introduced to satisfy the Einstein 
equations, whose form does not play a role in the following, and 



2fi(c) 



(52) 



The relation between the 4 and 10-dimensional Planck mass is as before: 



Mi = Ml 



^ = Mg,Vol(fia), 



(53) 



where Vol(f2e) is the volume of the compact manifold measured in units of M\o. 

We are interested in two limiting forms of the warped metric Eq. 1501 When there is no warping, or when the 
warping is everywhere very small compared to c, we recover the metric Eq. 1371 Alternatively, if there is strong 



warping e 



-4A (y) 



> c in a localized region of the Calabi-Yau, then upon the shift 



c — > c 



(54) 



and substituting with Eq. I51[ the metric becomes 



ds l 



e 2A °^ exp 



2 x(aQ 
3 M 4 



(g${*)dx»dx v + 2d j Bd li cdx fl dy^ + e~ 2Ao{ -^ {y)dy i dy j . (55) 



4 A short survey of the phenomenological benefits and shortcomings of various models can be found in [42l | . 

5 This analysis can also be relevant for the IIA vacua of |4.'I . which, given a restriction on the allowed fluxes, are related to IIB vacua by 
mirror symmetry. Other IIA compactifications leading to dS models also exist [4-li l45t | , see however [46| . More work on the difficulty of 
constructing type IIA dS vacua, including several no-go theorems, can be found in e.g. |47H49(I . 
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1. The scalar potential 

The four-dimensional effective field theory derived from type IIB compactifications on conformal Calabi-Yau man- 
ifolds has J\f = 1 supersymmetry. It is characterised by a Kahler potential K and a flux-induced Gukov-Vafa-Witten 
(GVW) superpotential W [6lj |. which determine the scalar potential V for the moduli [54ll55j. The effective theory has 
a no-scale structure, implying that the flux-induced potential is independent of all Kahler moduli except the volume 
modulus, which sets the overall scale of the potential. The GVW superpotential and no-scale structure survives the 
corrections due to warping [H3, EH, which can be encapsulated entirely in the shift of the field c as above. The 
dependence on the volume modulus x is therefore a simple scaling 

V = exp(-V6x/M 4 )V (<j>), (56) 

where <j> collectively denotes all the flux-fixed fields of the compactification. Thus x is a runaway direction in the tree 
level potential, exactly as in the toy model discussed in Section [TTJ We will see below that the volume modulus also 
couples to the branes in the theory, just as in the Einstein-Maxwell example, and will obstruct their stability if it is 
left unfixed. To fix this modulus, sub-leading quantum corrections must be added to the potential. 

Existing mechanisms for moduli stabilization, such as KKLT (5|| and BBCQ [62j , take advantage of non-perturbative 
g s or perturbative a' corrections that induce a potential for the Kahler moduli. In both cases, an AdS minimum can 
be produced. In the KKLT model, only the g s corrections are considered. The AdS minimum is supersymmetric and 
of depth 

a? A^e~ 2aCmin 

V K KLT{c m \n) ~ -z , (57) 

OC m i n 

where a and A are undetermined (but in principle calculable) constants appearing in the non-perturbative corrections 

to the superpotential W — Wo + Ae ac [56j]. 6 Since c m j n ~ — ^ log , and we require c m i n 3> 1, the construction 

only works for very small, but non-zero, vacuum expectation values of the tree level superpotential Wq. In the absence 
of other sources, the full potential approaches zero from below at large c. 

In the BBCQ model (a.k.a. LARGE volume compactifications) both g s and a' corrections are considered. To 
obtain a vacuum in this model, it is necessary to consider two different Kahler moduli, where the overall volume 
of the compactification is some linear combination of the two fields. For small values of the flux superpotential 
| Wq I , a KKLT minimum is reproduced, whereas for larger | Wq | a non-supersymmctric AdS minimum is obtained at 
exponentially large volume of depth [62( 

V BBCQ (Xmin) « -e K -\A s W \^-e- 3 Vlx/M^ (58) 

The Kahler potential K cs and flux superpotential Wo are in principle calculable from the fluxes and vacuum expec- 
tation values of the complex structure moduli and axio-dilaton. The constant A s arises from the non-perturbative 
corrections to the potential for the Kahler moduli. Again, the potential approaches zero from below at large volume. 

As already noted, both KKLT and BBCQ typically lead to AdS vacua, which must then be "uplifted" to Minkowski 
or dS vacua by adding extra sources to the theory. This can be achieved by space-filling anti-D3 branes localized in 
a warped throat of the Calabi-Yau [56[. Alternatively, one can use space-filling, partially wrapped Z?7-branes with 
world- volume gauge fluxes [63]. The scaling of the contribution of the anti-D3 branes with volume and warping is 
obtained from the action Eq. using the metric ansatz Eq. [50] in the appropriate limit: 



strong warping : S3 



rf 4 ^^, (59) 
weak warping : S3 — J d^x^/g — |, (60) 

where 173 is the brane tension, and we have neglected the effect of the 5-form flux sourced by the brane (roughly, this 
contributes a factor of 2 [56| ) . 

In the case of a partially wrap ped i?7-brane, the uplifting term in the four-dimensional action is proportional to 
the integrated world- volume flux [63l |64[, and scales similarly to the contribution from the anti-i)3 branes. Allowing 



The (possibly large) shift in c made above can be absorbed into the coefficient A as in [60 
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for the possibility of stacks of branes, these new ingredients add an uplifting term to the potential for the volume 
modulus of the form: 

D3 or D7 without warping : V upUft = Der^ x/Mi , (61) 

D3 or D7 with warping : V upHft = De , 

with the constant D dependent on the warping, volume of wrapped cycles, gauge dynamics, and number of branes. 
These sources change the asymptotics of the potential such that it approaches zero from above, and resemble the 
Uplifted type of potentials in the toy model of Sec. HU 

In order to obtain a vacuum after uplifting, Eq. [61] must be comparable in size to the other contributions to the 
potential at the location of the original AdS minimum (that is, Eq. [61] must roughly be of order the depth of the 
AdS minimum). This generally requires the coefficient D to be very small since the uplifting terms fall off much more 
gradually with increasing x ( or equivalently c) than the other terms in the potential. In the KKLT example, D is 
made small by locating a stack of anti-D3 branes at the bottom of a strongly warped throat. The anti-Z?3 brane 
action Eq. [SH] depends on the warping, which can lead to an exponentially small D. For BBCQ, it is possible to use 
D7 branes to uplift the AdS vacua [64], [65[. In this case there are a number of ways to obtain a small coefficient 
.D, including t he p resence of warping, features of the internal space, and contributions from the world volume gauge 
theory [gl,[63,[6g. It may also be possible for warped throats to exist in the BBCQ setup, opening up the possibility 
for anti-D3 uplifting, although a detailed analysis along the lines of Ref. [6(| is necessary to determine the viability 
of such a scenario. 



2. Domain Walls 



In string and M theory, there are (possibly wrapped) branes which are charged under the fluxes used to stabilize 
the moduli of the compactification. Moving across such branes, the fluxes are adjusted, leading to changes in the 
potential for various moduli. Depending on the compactification setup, this can lead to quite dramatic differences, 
such as changing the relative importance of various non-perturbative contributions to the potential for the Kahler 
moduli. For simplicity, and since it is a generic effect, we will here assume that the most relevant difference is a 
(small) change in vacuum energy. 

In Type IIB string theory, there are D5 and iV55-branes which can wrap internal 3-cycles of a Calabi-Yau com- 
pactification manifold and give rise to charged membranes in the four-dimensional effective theory. Again, there are 
processes by which these charged branes can be nucleated from a background configuration of the fluxes [67l - l69j . 7 
In a warped throat with D3 and anti-£>3 branes, the nucleation of wrapped NS5 branes can be described by the 
Kachru-Pearson-Verlinde (KPV) mechanism [72j (see also [73|-[75j]). In a more general setting, solitonic domain walls 
can be constructed between vacua for the flux-stabilized moduli fields [76j]. It was shown in [77} (see also e.g. 78 1 
and [l6|) that the tension of these solitonic domain walls have the correct scaling to describe a wrapped five-brane. 
In all of these scenarios, the four-dimensional vacuum energy can be different on either side of the brane. 

As in the Einstein-Maxwell case, we are mainly interested in the effect of the tension term on the dynamics of the 
volume modulus. The appropriate source action for an A55-brane wrapped on a 3-cycle of period IT3 (we define the 
periods in terms of the metric g^) in the four-dimensional Einstein frame is 



J sources 



4 V M £, 



n. 



s3/2 



v^yexp 



^voi(r2 6 

The tension of the five-branes in the original 10 dimensional action is [7 



&NS5 



(27T) 



3 x(0 
2 Ma 



(62) 



(63) 



For D5-branes, the tension is obtained by replacing g\^ 2 by gj 1 in the above formula. In all cases, we require that 
the proper volume of the 3-cycles e _3j4 °n3 3> M{q 3 . In a strongly warped region, we can have n3 <C 1 while still 
satisfying this constraint. 



Some caution is needed here, since it has been argued that brane nucleation is an intrinsically stringy process, lying outside of the 
supergravity approximation |7dl. |71| . 
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From Eq. [62l we see that the tension of the dimensionally reduced brane couples to the volume modulus. Thus, 
just as in the toy model of Sec. [TTJ x picks up a spatial derivative at the location of the brane. Assuming that x has 
a fixed value x* at the wall, the change in the spatial derivative of x is 



n. 



V Aff ; ^Vol^g) 3 / 2 , 





/3 X*~ 


j exp 


V 2M 4 



(64) 



As in Section IIII A[ we check if the domain walls are stable by comparing this 'kick' velocity to the rolling velocity 
at the location of the brane. We first investigate the stability of domain walls before the potential is uplifted, i.e. for 
domain walls separating the supersymmctric AdS vacua of the KKLT model or the non-supersymmetric AdS vacua 
of the BBCQ model. For both models, the potentials are physical realizations of the toy model 'AdS' potential of 
Eq. [1] and static domain walls can always be found. In the KKLT model, these planar domain walls preserve some 
supersymmetry, and correspond to the well-known BPS solutions discussed in e.g. [l6j]. In the BBCQ model, there 
is no supersymmetry guaranteeing the stability of the solutions, but nevertheless, by the arguments presented in 
Sec. IIII A1 we find that the domain walls are at least stable against decompactification. 

Proceeding to domain walls in the Uplifted potentials, the relevant test for a solution is again to compare the 
rolling velocity to the local maximum of the potential, and (in the case of planar walls) the asymptotic region x oo. 
Starting with the latter test, we recall that stable domain wall solutions are always possible if the asymptotic falloff 
of the kick velocity (Eq. is faster than that of the rolling velocity. If energy is conserved, and the minimum has 
zero energy, then the rolling velocity is given by iV^j/tl for both the KKLT and the BBCQ models. Consequently, 
for either the anti-D3- or Z?7-branes, the falloff of the rolling velocity is too fast to allow for stable planar domain 
walls. 

The more appropriate test is therefore to compare the kick to the rolling velocity at the location of the potential 
maximum, which will uniquely determine if walls are stable. The local maxima of the KKLT and the BBCQ potentials 
are different, and we analyze the models separately. In both cases, we approximate the height of the potential barrier 
with the depth of the AdS minimum, since this will approximately set the necessary size of the positive uplifting term. 

The stability of domain walls in the KKLT model was investigated in Ref. 0] . Using Eq. 021 we reproduce the 
condition for stability 

(Ax') „*„ m B 3 3 e- :<i (65) 



2v/2|U max | Mg, Vol(ft 6 ) 3/2 V2aA c 

For the KKLT construction to be valid, we must have c m ; n 3> 1. In fact, there is a relatively narrow window in which 
all of the approximations made in this construction are satisfied, estimated by Ref. 74] to be roughly 10 3 < c m i n < 10 5 . 
At such large volume, to satisfy the criterion for stability, some of the other parameters in Eq. l65l must be exponentially 
small. From Eq. 1631 the brane tension decreases at small string coupling, but obtaining an exponentially small string 
coupling would seem unnaturally fine-tuned (and perhaps in conflict with phenomenological constraints on the string 
scale). The constant a, if it arises from gaugino condensation on D7 branes where a = 2~k/N]-> 7, c annot be arbitrarily 
small since Ndi cannot be large without large backreaction. In some cases, A might be large [80jj , but the only truly 
viable candidate is II3, which in a region of strong warping scales like II3 ~ e 3A °. 

Therefore, the stability of branes in this situation relies crucially on the period of the cycle that they wrap. If 
the 3-cycle is in the same throat as the anti--D3 branes used to uplift the potential for the volume modulus, we can 
estimate the period. From Eq. [55] and [HIl we see that the uplifting term scales like e 4A °. Comparing with the depth 
of the AdS minimum Eq. 1571 we must have 

e 4Ao occ min e- 2ac ™. (66) 
Therefore, the period of the three-cycle scales like 

n 3 cxc^ 4 ne - 3ac ™ /2 , (67) 

which, inserting into Eq. 1651 is sufficiently small to guarantee stable domain walls at large volume. In conclusion, 
stable domain walls can be constructed only if the cycle wrapped by the brane is sufficiently warped. This result 
was already suggested in where a similar result for solitonic domain walls was obtained. It is also interesting to 
compare our result to the near-BPS domain walls in Ref. (l6| . where supersymmetry-breaking terms have to be small. 
This is guaranteed by a large warping in the KKLT model, which is the same condition needed for the domain wall 
to be stable against decompactification in our model. 

Similarly, the condition for stability in the BBCQ model is 



(Ax') <J NS5 n 3 M 4 1/2 e- K "/ 2 e^«^? 



« M« Vol(tt 6 ) 3/2 \A S W \^ X )H 



< 1. (68) 
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Here, the volume is exponentially large, e.g. in the example of [62j we have x-min ~ 20 (yielding a volume of about 
10 10 in units of Mio). Thus, to obtain stable domain walls, a significant fine-tuning of the parameters in Eq. IB51 is 
needed. Just as for the KKLT example, most of these parameters are constrained. The expectation value of the 
superpotential, \Wq\ is determined by the fluxes. \Wq\ need not be small in BBCQ compactifications, but can in 
general not be arbitrarily large, since that would require large fluxes which might violate the tadpole conditions of 
the compactification. Therefore, we will assume that |Wo| is generically of order one, keeping in mind that larger 
values can be allowed for specific models. The complex structure Kahler potential K cs is determined by the unwarped 
periods, and is also of order one. As above, we find that the tuneable parameters are A s and II3, and to make II3 
small we need to put the corresponding cycle in a strongly warped region. Thus, even though large warping may not 
be necessary to construct vacua in the BBCQ model, it could be used to tune the period of the cycle wrapped by the 
five-brane. This seems necessary for stable domain walls. 

Before we end the discussion on the stability of domain walls in the BBCQ model, there are two things that we 
need to investigate. First, the a' correction to the Kahler potential couples the volume modulus to the axio-dilaton, 
which is a flux-fixed field. To disentangle the equations of motion for the two fields we assume that the axio-dilaton 
is fixed during the evolution of the volume modulus. This is motivated by the hierarchically small contribution of the 
a' corrections compared to the flux-induced potential. We leave a more careful study of these effects for future work. 

Second, in order to construct a vacuum at exponentially large volume, it is crucial that the compactification 
manifold has at least two Kahler moduli [62]. In the analysis above, we only considered the dynamics of the volume 
modulus, since this is what couples to the flux-fixed fields and the branes. However, the BBCQ potential, and also 
the non-trivial metric on the Kahler moduli space, implies that there are couplings between the volume modulus and 
the other Kahler moduli. We must therefore check, by looking at the equations of motion, that it is consistent to 
keep the other Kahler moduli fixed while the volume modulus is rolling. Performing this analysis, we find that this is 
possible, as long as the volume is exponentially large during the full evolution of the volume modulus. We therefore 
restrict our analysis to this regime, and hope to come back to the possible interplay between the Kahler moduli in 
future work. 



VII. CONCLUSIONS AND OUTLOOK 



The dynamics of the volume modulus can have important implications for the picture of eternal inflation in the 
string theory landscape. In this paper, we have shown that only a rather restricted set of configurations of the internal 
manifold and sources can be connected across stable planar or spherical domain walls. 

The dilatonic coupling between the volume modulus and domain walls separating different four-dimensional vacua 
can completely overwhelm the stabilizing potential in the neighborhood of the domain wall. Using analytic methods, 
we have shown that non-singular static planar domain wall and SO(3,l) bubble solutions exist only when the energy 
scale associated with the tension of the domain wall is smaller than the energy scale associated with the potential 
maximum for the volume modulus. Rephrasing this condition, solutions do not exist when the tension of domain 
walls separating compactified vacua is larger than the tension of domain walls interpolating between the compactified 
and decompactified phases, indicating that it is energetically favorable for such domain walls to decompactify. 

Using a set of numerical simulations, we find that this is indeed the case. When non-singular solutions do not exist, 
domain wall and bubble solutions seed a growing region of the decompactified phase. In the case of bubble solutions, 
when the true vacuum has zero or positive energy, the entire bubble interior goes over to the decompactified phase. 
This completely removes the spacetime region containing the true vacuum. Regardless of the transition mechanism 
that produces the original region of true vacuum, the true vacuum phase will not exist at late times, and therefore 
cannot be populated by any such transition. 

In a spacetime without a cosmological constant, the presence of an unstable domain wall separating different four- 
dimensional vacua causes all of space to decompactify. If there is a cosmological constant, roughly a comoving Hubble 
volume of the background space will be converted to the decompactified phase at late times. In neither case does the 
spacetime evolve to a solution where different four-dimensional vacua coexist. This conclusion is independent of the 
creation mechanism, be it the formation of domain walls via the Kibble mechanism, or the nucleation of bubbles of 
true vacuum during eternal inflation. 

On the other hand, our analysis shows that domain walls separating AdS vacua are not troubled by this instability. 
Although we have not used it in our analysis, it is interesting to note that such domain walls can be BPS, in which 
case their stability is guaranteed by supersymmetry. Our analysis is sufficient to show that such domain walls are not 
unstable to decompactification, but there may exist instabilities in other moduli fields if they couple to the domain 
walls. 

The dilatonic coupling of the volume modulus can also play an important role in determining the outcome of a 
collision between two bubbles. Such collisions might leave detectable experimental signatures of eternal inflation, but 
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only if our observed cosmology can exist to their future. For stable bubble walls, the energy released in a collision has 
the ability to seed a region where space decompactifies. We have performed a number of numerical simulations, which 
indicate that this process is in fact quite efficient. Depending on the kinematics and properties of the potential, the 
decompactified region can accelerate into the bubble interior, rendering such collisions unobservable. Further analysis 
of these preliminary results will undoubtedly be an important element in determining the observability of bubble 
collisions arising in context of eternal inflation driven by the string theory landscape. 

We have investigated compactifications of Einstein-Maxwell and type IIB string theory, which both give rise to 
landscapes of vacua, domain walls, and possibly eternal inflation. Fundamental branes act as domain walls separating 
four-dimensional vacua in these theories, and possess a dilatonic coupling to the volume modulus. We find that 
domain walls between dS or Minkowski vacua can be unstable to dccompactification, as described above. 

The dilatonic coupling becomes especially problematic when the potential fixing the volume modulus is small 
compared to the tension of the walls. This is a common feature for both KKLT and BBCQ compactifications of type 
IIB string theory since the volume modulus is unfixed at leading order. In this case, we find that stable domain walls 
and bubbles arise only when the brane wraps a strongly warped cycle, leading to a reduced tension. The instabilities 
provide further motivation for finding dS vacua in compactifications where all geometric moduli are fixed at the same 
scale, as is the case in e.g. generic type IIA compactifications. These compactifications bear some resemblance to the 
Einstein-Maxwell model, where we find that stability depends on the dimensionality of the cycles that a fundamental 
brane wraps in the internal manifold. 

The existence of unstable domain walls separating different four-dimensional vacua leads us to conclude that only 
a restricted set of configurations of the internal manifold and sources can be populated by eternal inflation. However, 
our results do not exclude the existence of solutions where vacua of different dimensionality, separated by domain 
walls, coexist. Such a solution could contain four-dimensional vacua, separated by intervening regions where a different 
number of dimensions is kept large, as suggested in Ref. [l5j . Further exploration of these and other ideas will be 
necessary to determine the full dynamical role of extra dimensions in populating landscapes of vacua. 
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